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Abstract 

Solving the stochastic Landau-Lifshitz-Gilbert equation numerically, we investigate the effect 
of the potential landscape on the attempt frequency of magnetization in nanomagnets with the 
thin-film geometry. Numerical estimates of the attempt frequency are analyzed in comparison with 
theoretical predictions from the Fokker-Planck equation for the Neel-Brown model. It is found that 
for a nanomagnet with the thin-film geometry, theoretically predicted values for the universal case 
are in excellent agreement with numerical estimates. 
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I. INTRODUCTION 



Much effort has been expended in fabricating deep sub-micron patterned magnets with 
the thin-film geometry. From the scientific point of view, such a small magnet is a good 
model system to study basic magnetism via a direct comparison with an idealized theo- 
retical prediction. From the application point of view, a steady progress of the patterning 
technology for fabricating a smaller cell has led to magnetic devices such as spin- valve read 
sensors for the hard disk drive and magnetic random access memories (MRAMs) utilizing 
the spin-transfer torque (STT) [l, 2| into a higher density, i.e. a smaller magnetic volume. 

Thermal agitation of a magnetization becomes more and more important as the magnetic 
volume of a unit cell decreases. In the spin- valve read sensor, the so-called "mag- noise" is 
a manifestation of the thermally excited ferromagnetic resonance in the sensor stack 4|. 
In the STT-MRAM, thermal agitation hinders a continuous miniaturization of the device 
because it can cause spontaneous changes of magnetization direction from one stable state 
to another. 

Thermal relaxation time is a statistical time-scale for which a magnetization escapes 
from an initial local minimum state over an energy barrier. The thermal relaxation time r 
of a magnetization is described by the Neel-Brown model in the high energy barrier 

asymptote, r = /q^ exp[UB/kBT] where /o is the attempt frequency, Ub is the energy barrier 
measuring the difference between a local minimum and a saddle point, ks is the Boltzmann 
constant, and T is the temperature in Kelvin. 

Experimental studies on the thermal relaxation of magnetization generally assume a 
constant attempt frequency 3, Is]. However, Brown showed theoretically that the attempt 
frequency is not constant but depends on many parameters such as the damping constant 
and the magnetic properties Followed by Brown's initial work j^, theoretical formulae 
2f the attempt frequency fo. dffie.ent potential .y,nn,et.y we.e proposed Q 
15|,ll6|,ll7|. 
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Accurate theoretical formulae of the attempt frequency are necessary for modeling ex- 
periments and predicting quantitatively the superparamagnetic limit for device applications. 
However, it is not easy to experimentally verify the theoretical formulae because i) an exper- 
imentally measurable quantity such as the switching field is mostly governed by the energy 
barrier, not by the attempt frequency, and ii) the damping constant, a key factor affecting 
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the attempt frequency, of a small magnet is not definite in general |18l . Il9 |. 

In this work, by means of a numerical study based on the stochastic Landau-Lifshitz- 
Gilbert (LLG) equation [sl, Q, we investigate the validity of the proposed theoretical for- 
mulae. It is found that for a nanomagnet with the thin-film geometry, theoretically predicted 
values for the universal case are in excellent agreement with numerical estimates whereas 
theoretical values for the intermediate-to-high damping limit and the very low damping limit 
fail to reproduce numerical ones in practically meaningful ranges of the damping constant. 

This paper is organized as follows. After introducing the proposed theoretical formulae 
(Sec. II) and numerical model used in this work (Sec. Ill), we show in Sec. IV the effect of 
the shape anisotropy, i.e. potential landscape, on the attempt frequency for various damping 
constants and discuss about validity of the theoretical formulae by comparing with numerical 
estimates. In Sec V, we summarize this work. 

II. THEORETICAL FORMULAE OF THE ATTEMPT FREQUENCY 

The magnetic potential of a single-domain particle with uniaxial symmetry in the 
presence of a static external longitudinal field H is given by 

U = Ku{l - al) - MsHai + ^{Niaf + N2al + N^a^M^ 
= Kjj ( 1 — sin^ 'd cos^ if) — MsH sin cos ip 

+ ^(A^i sin^ § cos^ ip + N2 sin^ § sin^ <^ + cos^ ^)^f , (1) 

where Ku is the uniaxial anisotropy, Ms is the saturation magnetization, H is the external 
field applied along the magnetic easy axis, d and ip are the polar angle of magnetization 
vector and the azimuthal angle of magnetization vector, respectively, (i = 1, 2, 3) is the 
direction cosines of magnetization vector, and Ni is the demagnetization factor along axis. 

In an axially symmetric potential {U{i),ip) = U{i})), Brown 
frequency is not a constant, but a complex function as 



&\ showed the attempt 



where 7 is the gyromagnetic ratio, a is the damping constant, is the effective anisotropy 
field, V is the magnetic volume, h is H/Hk, and H is the external field applied along the 
magnetic easy axis. 




FIG. 1: (Color online) Magnetic potential surface of a single domain particle with non-axial sym- 
metry, ai is the magnetic easy axis and the external field H is applied along the easy axis. The 
magnetic energy U in the Eq. (1) has two equivalent saddle points and two minima points; local 
minimum and global minimum. In thermally activated switching, the magnetization changes from 
local minimum to global minimum, passing through saddle point. 



4 



Dependence of the attempt frequency on the damping constant for a non-axially symmet- 



ric potential was first theoretically predicted for two limiting cases; i 
damping (IHD) case [lOl , and ii) very low damping (VLD) case 



versal theoretical equation 
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intermediate-to-high 



12 1 . Later, the uni- 



171 ] which is valid for all values of the damping 



constant was derived by extending the Meshkov-Mel'ikov depopulation factor to the mag- 
netic case 2l|]. In this work, numerical estimates of the attempt frequencies are compared 



with the theories for the two limiting cases, and the universal one. 

n the IHD limit, the attempt frequency for a non-axially symmetric potential is given 
as lol |. 



Jo 



IHD 



ja^CmiCm2 (4 - Ci) + ,J (ci + c'a)^ + 4cic'2/a2 



(3) 



where Cmi and Cm2 are the coefficients in the expansion of magnetic potential U about a local 
energy minimum for the initial magnetic state, U = f/m + l/2(cmiQ;f + Cm2tt2) + ' ' ' J Ci and c'2 
are the coefficients in the expansion about the saddle point, U = f/g + l/2(cia^ — C2al) + ■ ■ •, 
respectively. 



Klik and Gunther 



ll[ |. and Coffey 12| derived a theoretical formalism of the attempt 



frequency for a non-axially symmetric potential in the VLD limit as, 

cVLD 



rv 

Jo 



■5, 



1'nMs{\ + a 

where S is the dimensionless action variable at the saddle point potential Us defined as 



(4) 



5" 



V 

ksT Ju{^,<fi)=Us 



;i -cos^ ^) 



d 



dcos-t} 



1 



d 



[1 — cos^ ^) dip 



U{^,^)dcos^ . (5) 



For the universal case, the attempt frequency is given by 



Q, Q, Q 
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17|, 



IHD 



fo = AiaS)fi 

where S is given by Eq. ([5]). 

A{aS) is a factor which interpolates between the VLD and IHD limits, and given by 



(6) 



A{aS) = exp 



1 l n[l - exp{-ag(A^ + 1/4)}] 
A2 + 1/4 



TC Jo 



dX 



(7) 
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FIG. 2: (Color online) Probability of not switching (l — Psw) versus magnetization switching time 
at a=0.03 in the sample A {{I x w x d = 21 x 20 x 20nm^). 

III. NUMERICAL MODEL 

We performed macrospin calculations by means of the stochastic LLG equation, 

^ = -7M X H.. + (8) 

where M is the magnetization vector, and is the effective magnetic field including the 
external, the magnetostatic, the thermal fluctuation. To estimate the thermal relaxation 
time r, we used a macrospin model with Us/kBT of about 10 because of excessive compu- 
tation time. Probability of switching Psw of thermally-activated switching was estimated 
by counting the number of successful switching out of 500 switching events. The attempt 
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frequency is obtained by fitting numerical results of Psw as a function of the time using the 
Arrhenius-Neel decay of the probability of switching, Psw = 1 — exp[— /ot exp(— f/^/ZcBT)] 
as shown in Fig. [21 



IV. EFFECT OF SHAPE ANISOTROPY ON THE ATTEMPT FREQUENCY 

We have calculated attempt frequencies of various sized nanomagnets; sample A {I xw x 
(i = 21 X 20 X 20nm^), sample B (25 x 21 x IQnm^), sample C (40 x 30 x 7nm^), and sample D 
(100 X 28 X Snm^) where / (length), w (width), and d (thickness) are the sample dimensions 
along X, y, and z-axis, respectively and thus the x-axis is the easy axis. For all four samples, 
constant values of volume V{= 8400nm^), effective in-plane anisotropy Hk{= 875. 40e), and 
external field H (= -540 Oe) were used to exclude their effects on the attempt frequency. 
The thermal stability factor Us/kBT was 10.425, a good number for the high energy barrier 
approximation. The magnetic potential U of the sample A is axially symmetric since w = d, 
whereas U of other samples are non-axially symmetric since w ^ d. 

The numerical results of the attempt frequency for the nano magnet with an axially sym- 
metric potential (sample A) are shown in Fig. [3t^a). To our knowledge, the Eq. ([2]) was 
tested once by adopting the same way used in this work [22], and it was reported that the 
theoretical value of attempt frequency is different from the numerically estimated value by 
an order of magnitude. This inconsistency may have prevented further numerical studies 
on the attempt frequency. However, we found excellent agreement between the Eq. (|2]) and 
numerically estimated values (see Fig. 3(a)). The difference between the result in the Ref. 



22l | and ours originates from the sign of h in the Eq. (j2]). The h should be negative since 



the magnetization is initially in a shallower local energy minimum, whereas it was assumed 



to be positive in the Ref. 221]. The excellent agreement verifies the validity of our numerical 
approach to estimate the attempt frequency in this work. 

Fig. [3t^b), (c), and (d) show the dependence of the attempt frequency on the damping 
constant for the sample B, C, and D, respectively. Two features are worth mentioning. First, 
the attempt frequency increases with increasing w/d. For instance, the attempt frequency of 
sample D is an order of magnitude higher than that of sample A in the wide range of damping 
constant (Fig. [3]). Second, when the potential landscape is non-axially symmetric, there 
are two regimes of the damping constant where the attempt frequency shows an explicitly 
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FIG. 3: (Color online)(a), (b), (c), and (d) show the attempt frequency as a function of the damping 
constant a for the sample A, B, C, and D, respectively. Solid lines are theoretically predicted values 
and symbols are numerical results. Dimensions of samples and parameters are given in the text 
{Ms = 800emu/cm^). 

different dependence on the damping constant. At low damping {a < 0.03), the attempt 
frequency increases with the damping constant whereas at high damping {a > 0.03), it 
shghtly decreases. 

Considering the increase of attempt frequency with w/d, it should be noted that both 
Eqs. (13]) and (jl]) contain ^ycmlCm2 which is an averaged curvature of potential at the local 
minimum. Other terms do not vary much with w/d. The Cmi is given by HkMs{1 + h) 
which is a constant for all four samples since both effective in-plane anisotropy field Hk and 
external field H are assumed to be constants. The Cm2 is given by HkMs{1 + h+ 
where N2 (N^) is the demagnetization factor along the in-plane hard (out-of-plane hard) 
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axis. The Cm2 significantly varies with the sample shape since is dominant. 

Therefore, an important parameter to determine the attempt frequency is the coefficient Cm2 
which measures the curvature of potential along the direction cosine a2 and is related to 
the out-of-plane demagnetization effect. A larger curvature of potential at a local minimum 
results in a higher attempt frequency. It is because the magnetization moving away from a 
local minimum due to a thermal random force experiences an instantaneous restoring force 
proportional to the curvature. The curvature Cm2 becomes smaller and smaller as the aspect 
ratio of sample w /d approaches the unity. Among the tested samples, the sample D provides 
the largest Cm2 and thus, the highest attempt frequency. 

In order to understand dependence of the attempt frequency on the damping constant, 
we compare numerical results with the theoretical formulae (Eq. ([3]) (IHD) , (j4j) (VLD) 
and (E]) (Universal)) from the Fokker-Planck equation for the Neel-Brown model. In the 
whole range of damping constant , the numerical results are in good agreements with the 
Eq. ([6]) multiplied by factor 2. The equations were derived for the escape of magnetization 
over only one shallower barrier assuming different barrier height between in-plane clockwise 
switching and counter clockwise one. In our case, the two energy barriers are identical since 
no symmetry breaking exists, validating the multiplication by factor 2. 

The theoretical values obtained from the Eq. ([3]) (IHD) partially coincide with the 
numerical results in high damping regime (a > 0.04), whereas the Eq. (jl]) (VLD) predicts 
much higher attempt frequencies than the numerical results in the tested range of damping 
constant (0.005 < a < 0.1). 

In the VLD, the escape rate is evaluated from the energy loss per cycle of a particle on 
the escape rate trajectory 15|, llTI]. The assumption made in deriving the Eq. (jlj), replacing 
the energy loss per cycle of the almost periodic motion at the barrier energy by the barrier 
height, is necessarily crude and only applies when the damping constant is less than about 
0.001. The failure of Eq. (jl]) to estimate the attempt frequency is also found in the Ref. [l4 | 
where comparisons among the IHD escape rates, the VLD escape rates, the universal solution 
based on the Meshkov-Mel'inkov depopulation factor, and the exact escape rate based on 
the continued fraction solution of the Fokker-Planck for the lowest eigenvalue were made. 
In the Ref. [l^, it is shown that the VLD asymptote begins to fail for the damping constant 
of the order of 10~^ even if the action on the escape trajectory is evaluated exactly whereas 
the universal solution provides a reasonably accurate approximation throughout the whole 
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range of damping. 

Therefore, it is obvious that the universal escape rate (Eq. ([6])) provides an accurate 
description of the behaviour of the exact escape rate provided that the barrier height is 
sufficient to allow one to define an escape rate. Furthermore, since the damping constant 



m a 
0.03 



Typical nanomagnet with the thin film geometry is in the range between 0.005 to 
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261] where the VLD and the IHD approximations show evidently wrong 



predictions for the attempt frequency, the Eq. ([6]) should be used to design experiments and 
to interpret experimental results performed at non-zero temperatures. 



V. SUMMARY 



In a nanomagnet with the thin-film geometry, the demagnetization energy along the 
magnetic hard axis is a a main factor affecting the attempt frequency. Comparing numerical 
estimates of the attempt frequency of magnetization with the theoretically predicted values, 
we verify the validity of the theoretical formula of the attempt frequency for the universal 
case. However, the theoretical formulae in the low damping limit and the intermediate- 
high damping limit fail to reproduce numerical values for the typical range of the damping 
constant. Therefore, the attempt frequency obtained from the theoretical equation for the 
universal case should be used to design experiments and to interpret experimental results 
performed at non-zero temperatures. 
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